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The causal dynamical triangulations (CDT) program has for the first 
time allowed for path-integral computation of correlation functions in full 
general relativity without symmetry reductions and taking into account 
Lorentzian signature. One of the most exciting recent results in CDT is 
the strong agreement of these computations with (minisuperspace) path in- 
tegral calculations in quantum cosmology. Herein I will describe my current 
project to compute minisuperspace (Friedman-Robertson- Walker) path in- 
tegrals with a range of different measures corresponding to various factor 
orderings of the Friedman-Robertson- Walker Hamiltonian. The aim is to 
compare with CDT results and ask whether CDT can shed light on factor- 
ordering ambiguities in quantum cosmology models. 

PACS numbers: 98.80.Qc, 04.60. Gw, 04.60.Nc 

1. Introduction 

In every approach to quantum gravity, persistent challenges from dif- 
feomorphism invariance rear their heads, with forms varying according to 
method of quantization. One particularly unregenerate challenge comes 
from the fact that the kinetic term in the Hamiltonian of a diffeomorphism- 
invariant theory contains products of position and momentum variables. In 
a canonical approach, this leads to factor ordering ambiguities because the 
kinetic term in the quantum Hamiltonian will contain products of noncom- 
muting operators, while in a path integral setting it results in a correspond- 
ing indeterminacy in the definition of the functional integral measure. 

Before even tackling the definition of a path integral measure for general 
relativity, however, we encounter a more immediate difficulty engendered by 
the sheer enormity of the set of paths to be integrated over. The space of 
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histories is through superspace, the set of all possible spatial geometries (3- 
metrics modulo diffeomorphisms) . In order to make sense of the integration, 
let alone perform practical computations, one is forced to restrict the class of 
geometries considered by imposing one or another simplifying assumption. 
Most common are two fundamentally different approaches to reducing the 
space of allowed geometries. 

Chronologically first came the method of reducing superspace to a "min- 
isuperspace" by restricting to geometries so symmetric that only finitely 
many degrees of freedom remain fH [2]- With the obvious convenience of 
reducing from a quantum field theory to an effectively quantum mechani- 
cal system, one can consider propagators and correlators for quantities of 
cosmological interest such as the scale factor. These advantages are al- 
loyed by the danger of constructing a quantum theory in which artifacts of 
the symmetry reduction overshadow genuine quantum gravitational physics; 
however, many issues to be faced in the quantization of full general relativity 
are preserved in these quantum cosmological models. As a mixed blessing, 
the twin problem of factor ordering and path integral measure persists. Be- 
cause of quantum cosmology's great potential power and utility, one would 
dearly like to know how far to trust it as a model problem for full quantum 
gravity. 

Another, newer approach to the reduction of superspace, the method 
of causal dynamical triangulations (CDT), imposes no symmetries but in- 
stead constructs a regularized path integral by restricting to piecewise flat 
geometries [3lUl[5]. Spatial slices are constructed of 3-simplices, and joined 
to adjacent spatial slices by edges which complete the whole into a 4-d sim- 
plicial manifold. These joining timelike edges are required to connect in 
a fashion which preserves the causal structure of the slicing (see e.g. 
so that although the path integral performed is Euclidean, only causally 
well-behaved histories are summed over. 

The Euclidean Einstein-Hilbert action is replaced in CDT by the Eu- 
clidean Regge action; the path integral by a sum over triangulations: 



where T is a causally well-behaved triangulation of spacetime (taken to be 
topologically I x S'^). The quantity Ct is equal to the order of automor- 
phism group of the triangulation T, and the 1 /Cq- factor in the summation 
is known as the discrete measure because it weights highly symmetric ge- 
ometries with a lower probability. In practice, the sum over triangulations 
is approximated using a Monte Carlo simulation which generates a set of 
independent histories. 
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Because the CDT simulation deals directly with (triangulated) geome- 
tries and computes discrete averages rather than continuum integrals, this 
choice of discrete measure to employ in the summation has a far more trans- 
parent meaning than is possible in quantum cosmological models. 

Recently, computations have become possible which allow comparison 
of results from CDT with those from the semiclassical or quadratic approxi- 
mations in the quantized Friedman-Robertson- Walker (FRW) cosmological 
model 0[7l[8]. The semiclassical state in FRW can be recovered from CDT 
simulations, and by computing the covariance matrix of the (discrete) spa- 
tial 3- volume, one can even observe CDT realizing quantum fluctuations in 
3-volume. To an impressive degree, these fluctuations match those given by 
the quadratic approximation of quantized FRW. 

It is in this second comparison, of the CDT covariance matrix with the 
FRW quadratic approximation for volume fluctuations, that I show how 
the signature of factor ordering can be detected. By performing an FRW 
quadratic approximation analogous to that in [7] and [8] but with varying 
choice of factor ordering and therefore varying path integral measure, I set 
up a range of possible path integral expressions with the goal of performing 
each and comparing to data from CDT to determine which quantum cos- 
mological path integral measure most closely approximates path integrals 
performed in a non-symmetry-reduced manner. In this way, CDT may 
be able to help us decide which path integral measure for FRW quantum 
cosmology most faithfully represents the full quantum gravitational path 
integral with all modes except the global scale factor integrated out. 

The plan of the paper is as follows. Section [2] summarizes the calcula- 
tion in |7j and [8j comparing CDT's 3-volume covariance matrix with the 
quadratic approximation for quantum fluctuations in the FRW model, and 
points out the stage at which path-integral measure dependence enters in. 
Section [3] adopts a 2-parameter family of factor orderings for FRW quantum 
cosmology first studied by Steigl and Hinterleitner in [9j , and uses a WKB- 
style approximation for the solution to the FRW Wheeler-DeWitt equation 
to estimate the computational effects at quadratic level of using different 
factor orderings. In Sectional I demonstrate that path integral expressions 
corresponding to all the Steigl-Hinterleitner factor orderings can be obtained 
in relation to the known path integral expression (due to DeWitt [lOJ and 
Parker [11]) for a particle in a potential on a curved background. Finally, 
Section [5] discusses the results and outlines directions of future research. 

2. CDT and the FRW model 

In the Monte Carlo simulations which provide the computational back- 
bone of CDT, measuring the covariance matrix for spatial volume is a con- 
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ceptually straightforward affair. Using the K independent histories gener- 
ated by the simulation, one keeps track of the spatial volume N^'^ {i) (the 
number of 3-simplices) for each history k, at each timestep i. Following the 
notation in [8], the expectation value for spatial volume is approximated by 
a simple average over the histories at hand: 

]V3(z)^(iV3«) = -ij]iVf (.). 

k 

The covariance matrix for spatial volume is now given by 

C {ij) = ^Y1 (^f (^) - ^3 (i)) (ivf (j) - Ns (j)) • (2.1) 

k 

To compare this covariance matrix with results from FRW quantum cos- 
mology, one must first write the Euclidean signature (Wick- and conformal- 
rotated) FRW Lagrangian in terms of the spatial volume rather than the 
scale factor: 

Denoting the quantum fluctuations in spatial volume by x{t) = V3(t) — 
V£\t), where V^\t) is the classical solution, the correlator for the vol- 
ume fluctuations can be approximated by the usual WKB expansion, to 
quadratic order: 

{x (t) X {t')) ^ ■ J Vx (s) X (t) X {t') e-27T If^' -is)M{s,s'Ms') 

= e ^ -M-^ {t,t') , (2.3) 

according to standard path integral computations. 

However the final equality in (j2.3p depends upon the assumption that 
the path integral is Gaussian; i.e. that taken together, the expression 
Vx (s) e~2h fl'^'^ x{s)M{s,s')x(s') -g functional Gaussian measure. Because 
the kinetic term in the Lagrangian <\2.2h involves both V3 and its time deriva- 
tive, the situation is in fact complicated by ambiguity in factor ordering 
and hence choice of path integral measure. Here we receive our first hint 
that perhaps even this simple quadratic approximation computation may 
be affected by the factor ordering problem. Far from being a misfortune, 
this added complexity in the quadratic approximation is our opportunity 
to probe the factor-ordering problem in quantum cosmology with currently 
accessible computational tools. 
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The goal of this paper is to revisit the computation (12. 3p . taking into 
account variation in path integral measure resulting from varying choice 
of factor ordering in the quantum FRW Hamiltonian. For ease of contact 
with other approaches to the factor ordering problem, the remainder of the 
paper considers the FRW Hamiltonian in terms of the scale factor rather 
than spatial volume, to which it is related simply by V3 = a^. 

I show that by ranging over the two-parameter family of factor order- 
ings for the quantum FRW Hamiltonian introduced in [9], the correlator 
computed in (j2.3p will vary in a well-defined and quantifiable fashion. Ul- 
timately the aim is to compare the computed correlators for each factor 
ordering with CDT data and determine a best fit. 

Before proceeding, a simple WKB approximation to the solution of the 
Wheeler-DeWitt equation for varying factor ordering gives us an assessment 
of the extent to which we may expect results at the quadratic level to depend 
on factor ordering. 

3. How much difference does choice of measure make? 

Classically, the Hamiltonian for the FRW universe is given by 



where as usual A is the cosmological constant. In [9], the form of the 
quantized Hamiltonian is allowed to range over a 2-parameter family of 
factor orderings, so that 



Two special cases in the family of Steigl-Hinterleitner factor orderings 
are worthy of note: first, by taking i = k, we obtain a 1-parameter family of 
symmetric factor orderings for which the Hamiltonian (13. 2p is self- adjoint 
with respect to the ordinary Lebesgue measure da. Alternatively, by con- 
sidering i = j = 1/2, k = 0, we can make the kinetic term of ()3.2p into a 
covariant Laplace-Beltrami derivative operator 




(3.1) 




(3.2) 



where i + j + k = 1 and V{a) = ^ ( 



a+ ^a^). This Hamiltonian is self- 





(3.3) 
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which for a 1-dimensional metric g{a) = (a) is given by 

-V2a,a-V2a,. (3.4) 



a 



This choice i = j = 1/2, which we can call the "Laplace-Beltrami" factor 
ordering, is the one we would choose in order to regard the Hamiltonian (j3.ip 
as effectively that of a particle in a potential V{a) on a curved background 
g{a) = (a). 

The quantized Hamiltonian (|3.2p yields the Wheeler-DeWitt equation 



V'(a) = 0, (3.5) 



h^a-^dl^p (a) + h^{-k + i- 1) a-'^dai^ (a) 



+ 



h^k i-i + 2) - a + 



3 

for which we can construct a WKB-style series solution of the form 

(a) = exp |-^| f; n>„ (a) . (3.6) 

71=0 

In standard path integral treatments, one uses arguments from the method 
of steepest descent to assert that in the expansion of a path integral about 
its classical solution 

1 Iff 5^s . ,2 1 fn . ,3 



/l>.e>:p{-i(5. + i// 



6x1 ' 3! 777 (5x3 



(3.7) 



terms of order (x — XdY higher in the exponent yield terms in a summa- 
tion expansion of order h and higher. Using this reasoning, the quadratic 
approximation for the path integral is equivalent to the limit h ^ and 
thus to the truncation of the WKB series solution of the Wheeler-DeWitt 
equation from the form (j3.6p to 

V9o(x)exp|-^|. (3.8) 

The arguments necessary to make such an equivalence of approximations 
precise would depend on the definition of the path integral measure. In our 
case, this is precisely what is in question, so truncating the WKB series 
solution from (j3.6p to (j3.8p must not be taken as anything more than a very 
rough heuristic guide to what we may find upon performing a quadratic 
approximation of the path integral. Still, we shall find that solving for (j3.6p 
up to ipo proves to be an illuminating exercise. 



ordering printed on October 23, 2009 



7 



Inserting ()3.6p into the Wheeler-DeWitt equation (13. 5p yields, to second 
order, 



V30 (a) = K 



A 2 



k — i 

a—, 



(3.9a) 
(3.9b) 



where coefficient K in ()3.9bp is an arbitrary constant. As expected, (|3.9ap 
is the solution to the imaginary-time Hamilton-Jacobi equation. Hence 
exp [—5(a)] is the semiclassical state and is independent of factor ordering 
since (|3.9ap does not involve i,j, or k. 

As a clue to whether the quadratic approximation may carry the signa- 
ture of factor ordering, the important feature of (|3.9bp is the factor a~2~, 
which indicates that at quadratic level we may hope to measure the devia- 
tion in factor ordering from the symmetric case i = k. Emboldened by this 
result on the canonical side, we turn to path integral computations. 



4. Path integrals for Steigl-Hinterleitner orderings 

To compare our results with those from CDT simulations, we must con- 
sider correlators in the path integral formalism, in particular correlators 
for functions of the FRW scale factor. As in CDT, proper time will be 
used throughout for the FRW model. However to define path integrals with 
various measures in terms of time-slicing, it is more natural to deal with 
propagators than correlators. This mismatch is circumvented by the simple 
observation that correlators can always be broken down into expressions 
involving propagators: 

{a",t"\fiait2))f{a{h))\a',t')ijk = 

POO 

/ a'2^da2 f {a2) {a" ,t"\a2,t2)ijk x 
Jo 

' roo 

/ a'^^daif {ai) {a2,t2\ai,ti)ijk{ai,ti\a' ,t')ijk , 
Jo 

(4.1) 

by using the appropriate resolution of identity. The subscript ijk labels 
the correlators and propagators as those corresponding to a given choice of 
Steigl-Hinterleitner factor ordering of the FRW Hamiltonian. 

Since correlators can be expressed in terms of propagators, the focus can 
shift to finding a path integral expression for the propagator Kijk (a, t; a' , t') 
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satisfying 



^ - ^a-'daa-^daa-'' - V{a) ) K.^k (a, t; a',t') =0, a ^ a' 



lim Kijk (a, t:a ,t 



5 (a — a 



(4.2) 



For the case of the Laplace-Beltrami factor ordering, the path integral 
expression for this propagator is known [10\ [TT] . It is the propagator for a 
particle moving in a potential V{a) on a curved background g{a) = (a): 



Klb {a" ,t";a' ,t') 





a 2 {t)a (t) 




f 




/C{a",t"\a',t'} 





1 



aa^ — V{a) 



= lim TTT- 

N 



N-l 



X exp 



-T 



n=l 



n=l 
1 

2e 



/ andUn 



(4.3) 



In the remainder of this section, I show that path integral expressions for 
all other Steigl-Hinterleitner factor-ordered FRW Hamiltonians can be ob- 
tained in relation to (j4.3p . 

For any ordering (i, k) in the 2-parameter family of Steigl-Hinterleitner 
orderings, a two-step process relates the (i, k) propagator to (j4.3p . First, the 
propagator for any ordering belonging to the line i + k = 1/2 can be related 
to ()4.3p by convolution with appropriately chosen Green's functions. Next, 
the propagator for an ordering belonging to any perpendicular line i — k = c 
(c = constant) can be related to a propagator whose kinetic term is ordered 
according to i — k = c, i + k = 1/2, and whose potential term differs from 
V{a) by the addition of a quantum potential proportional to h?. Using the 
convolution result described above for orderings on the line i + k = 1/2, 
this new propagator with added quantum potential can then be related to 
the propagator (j4.3p with a corresponding added quantum potential. These 
results are described in detail below. 



4-1. Orderings with i + k = 1/2 
For orderings lying along the line i + k = 1/2, the quantum FRW Hamil- 



tonian is 



H 



(4.4) 
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so we need a propagator satisfying 

- Kijk {a", t"; a, t') = 0, a" / a (4.5a) 

Ihn Kijk [a" ,t";a ,t') = (a')^^ ^^'^5 (a" — a) . (4.5b) 

Note that we can write the Schrodinger operator in (j4.5ap as 

^ H)=a''(i-HLB]a-^ (4.6) 



dt J \dt 

where Hlb is the quantum FRW Hamiltonian with Laplace-Beltrami factor 
ordering. Hence by convolving the propagator Klb for the Laplace-Beltrami 
Schrodinger operator with the Green's functions a^S{a — a') and a~^5{a — a') 
for the multiplication operators and , and converting to the correct 
measure a}/'^~'^^da^ we obtain the propagator Kij]^: 

K,^u {a",t";a',t') = {af'-'/' {af/' [{a'fKLB {a" , t" ; a' , t') (a')"'" 

= {a'a'fKLB{a",t";a',t'). (4.7) 



4-2. Orderings with i — k = c 

All orderings lying outside the line i + k = 1/2 are dealt with according 
to which perpendicular line i — k = c they occupy. First, rewrite the quan- 
tum Hamiltonian in terms of a new effective kinetic term and a quantum 
potential added to V{a): 



H 



V{a) 



2 
2 



[a-^dl - (1 - c) a-^da + k{2-c-k) a'^] + V{a) 



^^l/2(l/2+c)^^^-l/2^^^-l/2(l/2-c)- 



k{2 



' 4 12 



CM--C 



+ Y{a) 



(4.8) 



Notice that the effective kinetic term in the second to last line is ordered 
according to a different Steigl-Hinterleitner ordering having % = (1/2 -|- 
c)/2, / = 1/2, A;' = (1/2 - c)/2. This ordering satisfies i' /c' = 1/2, so 
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now the convolution procedure of Section [4.11 applies . with the new effective 
potential on the last line of (j4.8p instead of simply V{a). Thus we obtain 



Ki,k {a",t";a',t') = {a' a") 

^1 



1/4- 



k — i 
2 



V 



a2 



exp 



dt 



-aa 



C{a",t"\a',t'} 

y(a) + yp(/c,c) a"^ 



(t)a(t) 



(4.9) 



p{k,c) = k{2-c-k)-U^ 



--C 



The expression (|4.9p gives the path integral expression for the propagator 
corresponding to any Steigl-Hinterleiter ordering, as a well-defined deviation 
from the Laplace-Beltrami case. The deviation from the Laplace-Beltrami 
propagator enters in two ways: through the quantum potential measuring 
the ordering's distance from the line i + k = 1/2, and through the prefactor 
measuring the ordering's degree of asymmetry. 

Reversing the breakdown of correlators into propagators given by (|4.ip . 
we can use (j4.9p to express any Steigl-Hinterleitner correlator for a function 
/(a) of the scale factor in terms of the Laplace-Beltrami correlator for the 
same function: 



{a",t"\fia{t2))f{aih))\a',t%k 

' ' ^ (a", t"\f (a (t2)) / (a (ti)) \a' , t')LB,QP, 



(a"a'' 



(4.10) 



where the subscript LB, QP denotes the fact that the correlator on the right- 
hand side is taken with the Laplace-Beltrami factor ordering, and with the 
added quantum potential. 



5. Discussion 

As a confirmation of the WKB prediction in Section [3l notice that the 
prefactor 



{a"a'y 



in (j4.9p and (j4.10p measures the effect on propagators and correlators of the 
factor ordering's degree of asymmetry, in a manner closely resembling that 
of the ordering-dependence in the first term of the WKB fluctuation factor 
(Hlbl): 



<po (o) = K 



A 



k — i 

a 2 



(5.1) 
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The obvious next step is to compute the path integral in (14. 9p to quadratic 
order, expecting based on WKB analysis that ordering dependence will per- 
sist. To compute the quadratic approximation to the path integral, tech- 
niques such as those applied in [12] and references therein may prove useful. 

Having computed the path integrals in (j4.9p and (j4.10p to quadratic 
level, the goal is of course to compare with the CDT data for volume fluc- 
tuations described in Section [2l This will necessitate the conversion from 
results in terms of the scale factor to results in terms of spatial volume. 
However, this is a straightforward affair which could be approached either 
by computing correlators in terms of a? or by recreating the above deriva- 
tions for the Lagrangian (j2.2p . 

Once we have computed the FRW volume fluctuation correlator for each 
path integral measure [Da] in the 2-parameter family, we can compare 
against the CDT covariance matrix (j2.ip and determine the (i, k) which 
yields a best fit. The covariance matrix (j2.ip is computed for a fixed number 
A^4 of 4-simplices constituting the triangulations; thus in line with standard 
finite-size scaling techniques, we should compute (j2.ip for a sequence of 
increasing A^4, at each stage finding a best-fit (i, k)i\f^ until the fit stabilizes 
at some (z, k) = limAr^^oo(^, k) . The path integral computed with measure 
[Pa]^^ should then represent the continuum limit of CDT, with all modes 
integrated out except for the scale factor. 

An extension of the work presented here is to derive the path integral 
propagators for all Steigl-Hinterleitner orderings directly rather than in re- 
lation to the Laplace-Beltrami case. While the latter method (that followed 
in the current paper) may be more efficient for computing a comparison of 
correlators across operator orderings, a direct derivation of path integral ex- 
pressions for all Steigl-Hinterleitner orderings is likely to afford conceptual 
insight into the varying path integral measures. The aim is to construct 
such a derivation along the lines of that for the Laplace-Beltrami case using 
a quantum correction from the Weyl ordering prescription and short-time 
propagator [IB] . 

Further avenues for exploration are comparisons of the CDT data with 
computations in more sophisticated quantum cosmological models such as 
Bianchi IX [14j. Allowing the factor ordering of Bianchi IX to vary analo- 
gously to the Steigl-Hinterleitner orderings for FRW, it would be interesting 
to ask whether the same ordering for the kinetic term of the scale factor 
yields a best fit in both cases. More generally, consider a sequence of cosmo- 
logical models such that the degrees of freedom included in a given model 
are a subset of those in the next model (e.g., FRW and Bianchi IX). By 
comparing results from such a sequence of models with computations from 
a CDT simulation, one could ask whether the best fit factor ordering for 
a given mode (e.g. the scale factor) remains constant across cosmological 
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models of increasing complexity. In this way, it would be possible to ask 
what is the role of factor ordering in quantum cosmological models: does the 
best-fit ordering change at each stage, simply compensating for other modes 
which have been neglected/integrated out? Alternatively, does the order- 
ing remain stable with the inclusion of new modes, suggesting that it may 
genuinely reflect qualities of the factor ordering for full general relativity? 

To form an estimate of the insights likely to be gained from such a se- 
quence of cosmological models, one could study an analogous sequence of 
cosmologies in two spacetime dimensions, where CDT's continuum measure 
is analytically known (see e.g. [15] for coverage of 2d CDT). This line of 
inquiry may lead to hypotheses in an important open area: the nature of 
four-dimensional CDT's continuum measure. 
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